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Abstract 

We compute the first two moments {{sh — fno'^Y''^) of the final hadronic 
invariant mass in the inclusive decay B — > in the presence of a cut E"™™ 

on the charged lepton energy. These moments may be measured directly by 
experiments at the T(4S') using the neutrino reconstruction technique, which 
requires such a cut. Measurement of these moments will place constraints on 
the nonperturbative parameters A and Ai, which are relevant for extracting 
the quark masses m;, and rric, as well as the CKM angle Vcb- We include terms 
of order a^/3o and l/mf in the operator product expansion, and use the latter 
to estimate the theoretical uncertainty in the extraction of A and Ai. 
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I. INTRODUCTION 



Semileptonic B decays are a rich laboratory in which to study the bound state structure 
of the B meson. Of particular importance are the nonperturbative parameters A, Ai and A2, 
which arise in the heavy quark expansion for the meson mass ms to relative order 1/ml, 



\ _|_ 'iX 

ms = m, + A~ ^^—^ + 0(A|cdK) , (1-1) 



where 



Ai = l,B\K{%DfK\B) l2mB , 

A2 = (5| hJ^-a^^'G^.K \B)/QmB • (1.2) 



Here is the heavy quark field, and is the b quark pole mass. Since these quantities also 
parameterize the inclusive semileptonic and radiative decay rates of the B meson to order 
[^^,^,^,^, an accurate determination of them is essential for a reliable extraction of the 
CKM angle \ Vcb\ from inclusive semileptonic B decays. The parameter A2 is the expectation 
value of the leading operator which breaks heavy quark spin symmetry and therefore may 
be determined from the B-B* mass difference, yielding A2 ~ 0.12 GeV . However, A and Ai 
cannot be determined solely from mass measurements. 

There has been much recent interest in using inclusive observables other than the full 
semileptonic decay width to extract A and Ai. An analysis of the decay based on the heavy 
quark expansion yields an expression for the doubly differential decay rate dV /dq^d{v ■ q), 
where is the total momentum of the leptons and is the four-velocity of the B meson. 
Any observable which may be constructed from this doubly differential rate is sensitive to a 
linear combination of A and Ai. Actually, the simplest differential distribution to measure 
experimentally is the energy spectrum of the charged lepton, dT/dEi, which is somewhat 
more complicated because it depends on more than just q^. Study of this distribution has 
already yielded useful constraints on A and Ai [§|J3J^,0]- 

In two previous publications [|1C], we suggested that moments of the hadronic invariant 



mass in semileptonic B decay would also be interesting to study. If sh is the invariant 
mass of the hadrons produced in the semileptonic decay, and Wd = {uid + 3m£)*)/4: is the 
spin-averaged D meson mass, then positive moments of sh—^'d vanish in the parton model 
at tree level. Hence they are particularly sensitive to the power corrections proportional 
to A and Ai. However, it is difficult to measure these moments directly. Until recently, 
only fairly weak bounds on {{sh ~^'d)^) could be obtained, by combining information on 



various exclusive decay channels ||T^. With the introduction of the technique of neutrino 



reconstruction, this situation is changing [|ri[] . Soon it will be possible to measure sh directly 
and inclusively, by reconstructing the neutrino momentum pj^ and using it to find and 
V ■ q. The only complication is that this reconstruction requires a number of constraints on 
the phase space of the leptons, most importantly a lower cut on the charged lepton energy 
of about 1.5 GeV 0. 

While one might consider extrapolating the data to lower Ei and measuring the moments 
of Sh that way, it is clearly preferable to update the computation of the moments of sh to 
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include this lepton energy cut ab initio. It is the purpose of this note to present the results 
of such an analysis. There are few new theoretical issues which arise in this case, although 
the actual calculation is considerably more complicated than without the cut. We will refer 
the reader to our previous papers for a more complete exposition of the theory [jlO|, and 



only stress points where the new analysis differs in interesting ways from the old. We will 
also extend our earlier analyses by including the complete power corrections to the moments 
to 0{l/ml). Dimensional estimates of the size of these corrections will help us to estimate 
the error in the extraction of A and Ai due to 1/mf effects. They will turn out to be 
under control for the first moment, but very large for the second moment, compromising its 
usefulness for obtaining constraints on A and Ai. 



II. THE CALCULATION 

A. Power Corrections 

The analysis of the moments of sh begins with the doubly differential spectrum 
dV /dcpdiy ■ q), where = /ml and v ■ q = v ■ q/rrib. This distribution may be writ- 
ten as a product of a hadron tensor T^^ and a lepton tensor L^^ 

-^^J^ = r-(.,^)V(«..). (2.1) 

Each of these tensors has a decomposition in terms of scalar invariants, such as 

-le^'^^^v^qpL^ . (2.2) 

The Li are scalar functions of (f and v ■ q. Here we neglect the masses of the leptons. If we 
impose a cut Ei > xirtb on the charged lepton energy, then the effect is simply to modify 
the components Lj, 



^2 



(^4(f • qf - 3fv ■ q + 4[(f ■ qf - q^f^^ + 6xf 
-12x{v ■ qf + I2x^v -q-Sx^ 



247r[{v ■ g)2 - g2]3/2 



= 8n[{v ■ qf - q^f/^ ^'^'' - ■ g)(g^ - 4x.; ■ g + Ax') , 

The cut may be removed by setting it equal to the minimum charged lepton energy, x 

■ qf — q'^Y^'^, in which case Li = g^/Svr, L2 = = 0. Expressions for the analogous 
hadron tensor components Tj may be found in Refs. [0,^,^. The hadron tensor is independent 
of the cut X. 
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Since the hadron tensor T'^^ is computed with operator product expansion techniques 
which assume parton-hadron duahty, the calculation must be smeared by integration over at 
least one of or v ■ q before meaningful observables may be extracted . A nonzero cut 
X has the effect of restricting the available phase space for the leptons, which controls the 
range of integration in and v -q. If x is too large, the integration is too restricted, and the 
operator product expansion breaks down. This is known to happen, for example, near the 
lepton endpoint in charmless semileptonic B decays, for Ei above about 2.2 GeV [0,^,|TB|,|T3[ . 
However, our calculations of the coefficients of the 1/mf corrections indicate that the lepton 
energy cut of 1.5 GeV which is required for the neutrino reconstruction technique is not 
severe enough to cause such problems. 

The next step is to compute parton level moments of the form {E^{so — m^)"), where sq 
and Eo are the invariant mass and total energy of the strongly interacting partons produced 
in the semileptonic decay of the b quark. The energy is computed in the b rest frame. The 
partonic variables are related implicitly to the hadronic variable sh by 

Eq = 1- V ■ q, Sq = 1 ~ 2v ■ q + q"^ ; 

Sh = rn\ ~ "^rnBV ■ q + q^ ■ (2.4) 

Since these expressions involve both ms and rrtb, they must be inverted order by order in 
the heavy quark expansion using Eq. ( |1 . 1| ) . The final result will be moments of the form 
{{sh — ^d)™), for m = 1,2. In fact, we will see that only the first moment is really reliable, 
where the terms of order Aqcd and Aqqj-, in the heavy quark expansion are known. By 
contrast, the second moment starts only at order Aq^j-,, and hence is extremely sensitive to 
the large number of unknown parameters which arise at order Aq,-.^,. 



B. Radiative Corrections 



While it is possible to calculate the radiative corrections by calculating the 0{as) contri- 
butions to the Tj's themselves, it is much simpler to calculate directly the leading corrections 
to the parton model rate. The only subtlety in the calculation arises in determining the 
boundaries of phase space when the electron cut is imposed. Since the limits of integration 
of the electron energy depend on the parton level invariant mass Sq and lepton invariant 
mass squared g^, for given values of Sq and q"^ the electron cut may either lie below the lower 
limit of integration, in the region of integration, or above the upper limit of integration. Let 
/(■50)Q'^) be any smooth weighting function. Then the phase space integral is divided into 
kinematic regions, depending on the values of x, sq and g^. 



For X < 1(1 



rrir. 



dSn 



23:(1 — 3Q— 2a;) 



dq'^ds, 







+ / dso 

J(l-2x)2 JO 



l-2x 
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For X > |(1 — rric), 



.1_9t, 2x(1-so-2x) , 



(/(so,g )) = I dsQ dq f{so,q 



mt 



dq'^dsQ 



Here dVi/ dq^dsQ is the differential rate calculated with the electron energy cut imposed, and 
dV2/ dq^dsQ is the differential rate calculated with no cut, corresponding to x lying below 
the lower limit of integration for the electron energy. In the regions of phase space omitted 
from the expressions above, the cut lies above the upper limit of integration for the electron 
energy. 

The "BLM -enhanced" two-loop corrections are those which are proportional to 
af/^o, where (3q = 11 — Infj?) is the first term in the QCD beta function. These corrections 
dominate the two-loop corrections to many processes in QCD, and their effects in h decays 
have been discussed extensively in the literature |[T^. They are straightforward to calculate 



numerically using the techniques of Ref. [|T^], and no new subtleties are introduced into the 
calculation when an electron energy cut is added. Because of the renormalon ambiguity in 
its definition |]l8l, A is only defined order by order in perturbation theory. Since we are 



including the af/3o terms in our extraction of A, the resulting value is the "two-loop" A, and 
should only be compared with other extractions of A at the same order. 



III. RESULTS AND DISCUSSION 



The expansion of a moment of su takes the following general form, up to terms of relative 
order 



n 



+ 



'-'22 PO + <-^21 



niB 



A 



D. 



(m) 



A^ 



20 =2 



D 



(m) 



Ai 



B 



21 =2 

mi 



D 



(m) 



A, 



22 =2 



m 



B 



(3.1) 



+ 



A' 



(m) AAi 



'30 



m 



31 



D 



B 



m 



(m) 
32 



AA, 



(m) Pi 



B 



m 



33 



D 



B 



m 



B 



(m) P2 
34 — 3 



U ^B 

All the coefficients which appear are functions of the lepton energy cut i?™™. The parameters 
Pi and p2 are expectation values of local operators of dimension six which arise at order 1/ml 
in the heavy quark expansion. 



(3.2) 



where (if/ = 3 and dn = —1, respectively, for matrix elements between B and B* states. 
The % are related to nonlocal time ordered products of 1/ml terms in the operator product 
expansion with l/rrif, terms in the Lagrangian p4| . 
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hy i 



d^x / dt Ci{x)\B^*'>) +h.c. = 



oo 




{B^*^\K-a^^G^"' Ki f d^x f dt Cj{x)\B^*'>) +h.c. 

2 J J —oo 



Ti + dn^ 

rub 
% + dH% 
rub 



(3.3) 



The parameters pi and % are determined by nonperturbative QCD, and their values are not 
known; however, we compute their coefficients to ensure that none of them are anomalously 
large. We will use dimensional analysis to estimate pi and %, to obtain a rough estimate of 
the error in the extraction of A and Ai induced by 1/mf effects. 

The coefficients Dij are themselves functions of a^, although we will not compute radia- 
tive corrections to any of these coefficients. At present, Di^^ are known to order as only 
in the absence of an electron energy cut. While it would certainly be desirable to include 
terms of order asA/nib for general i?™™, the calculation is quite difficult and we have not 
attempted it here. This omission is particularly important for the second moment, since 
D^i'^ = 0. Without the cut, the term in question is numerically as large as the leading terms 
proportional to cj^"* and D^q (In the case of the ffist moment, by contrast, the radiative 
correction to D^^ is only a few percent |10|.) For this reason, as well as because of the large 
1/mf corrections, it is dangerous to use the second moment in a measurement of A and Ai. 
We will include our ignorance of these terms in the estimate of the theoretical error in our 
results below. 

Let y = Ef^^'^/rriB be the scaled lepton energy cut. When y = 0, we reproduce the known 
,0, 



results 
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The leading radiative correction to D'^i' is 0.099 Q;s(m6)/7r, and to is 0.038 Q;s(mf,)/7r 
For the preferred experimental cut of 1.5 GeV, for which y = 0.28, we find 
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The authors of Ref. M correct a numerical error in this term in Ref. |10|. 
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We present plots of the coefficients Cj and Di for arbitrary i?™™ in Figs. |l|-^ Note that 
the power corrections tend to blow up as the electron cutoff approaches its maximum value; 
fortunately, for a cut of 1.5 GeV the coefficients are not dramatically larger than without a 
cut. 

The effects of the l/mf corrections to the moments in the extraction of A and Ai are 
displayed in Fig. |. For the purpose of illustration, we ffist assume perfect experimental 
measurements of (sH — frlj^) = 0.30 GeV^ and {{sh — rnj^y) = 0.96 GeV"^, and then extract A 
and Ai using the theoretical expressions at order l/ml. We take as{mi,) = 0.22. These values 
for the ffist and second moments restrict A and Ai to lie on the solid and dashed curves, 
respectively, meeting at the point (A, Ai) = (0.31 GeV, —0.16 GeV^). The hypothesized data 
have been chosen so that this intersection coincides with the central values obtained from an 
existing analysis of the lepton energy spectrum [0. We then may estimate the theoretical 
uncertainty in our hypothetical result by following the approach of Ref. p|. By dimensional 
analysis, the parameters pi and % are all of order Aq(-;j3, and an estimate of their effect on 
the extraction of A and Ai is obtained by varying their magnitudes independently in the 
range — (0.5 GeV)^. Since the vacuum saturation approximation suggests that pi > 0, 
we take it to be positive, and we eliminate p2 by making use of the relation between p2, 
T2 and T4 and the D*-D and B*-B mass splittings presented in Ref. Finally, we vary 
the unknown coefficient of the asA/rrib term for both moments between half and twice its 
value with the cut removed. Varying the unknown parameters randomly in the allowed 
ranges, we find that the shaded ellipse shown in the figure and centered about the mean of 
the distribution contains 68% of the points. (Because the A^, AAi and AA2 terms bias the 
determination of A and Ai in a known way, the distribution is not centered about the point 
extracted from the theory at order 1/m^.). The region inside the ellipse gives a reasonable 
estimate of the theoretical error in the extraction of A and Ai due to higher order effects. 
Note that, as expected, the constraints from the second moment are strongly affected by 
higher order corrections, whereas the constraints from the first moment are quite tightly 
distributed about the leading result. Thus, as discussed earlier, only the hnear combination 
of A and Ai given by the ffist moment is significantly constrained. 

Constraints on A and Ai also have been obtained from moments of the lepton energy 
spectrum above 1.5 GeV 0]. To compare the theoretical errors in this approach to ours, we 
have performed an analysis analogous to that of the previous paragraph. This is similar to 
what was done in Ref. but we also include the terms proportional to a^Po The result is 
shown in Fig. ^. The size of the ellipse from the lepton energy analysis is slightly larger than, 
but comparable to, that which we obtained from the hadronic mass moments. Of course, the 
relative position of the ellipses is meaningless, since the hadronic mass moments have not 
yet been measured. Unfortunately, the two experiments effectively constrain the same linear 
combination of A and Ai, so that the measurements cannot be combined to determine both 
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parameters simultaneously. Instead, an observable sensitive to a different linear combination 
of A and Ai, such as the first moment in the photon spectrum in i? — > X^'^ [19|, will be 
required. In the meantime, consistency of the various allowed regions in the A-Ai plane will 
provide a powerful test of the heavy quark expansion for these decays. 



IV. CONCLUSIONS 



We have presented a calculation of the first two hadronic invariant mass moments in 
semileptonic B decay, in the presence of a moderate cut on the energy of the charged 
lepton. We included effects up to order alPo and 1/mf. These moments may be used to 
measure a linear combination of the HQET parameters A and Ai, with a theoretical accuracy 
which is comparable to, or slightly better than, the accuracy obtained from an analysis of 
the charged lepton energy spectrum. The consistency of the results obtained from these 
approaches will provide a test of the heavy quark expansion as applied to semileptonic B 
decays. To extract A and Ai simultaneously, it will be necessary to combine this analysis 
with that of a quantity sensitive to a different linear combination of the two parameters. 
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FIG. 1. Coefficients C^j' of the radiative corrections to the first and second moments, as 
functions of the electron energy cutoff in GeV. 
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FIG. 2. Coefficients of the power corrections to the first moment of the hadronic invariant 
mass spectrum, as functions of the electron energy cutoff in GeV. 
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FIG. 3. Coefficients D^j' of the power corrections to the second moment of the hadronic in- 
variant mass spectrum, as functions of the electron energy cutoff in GeV. 
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FIG. 4. Estimate of the theoretical uncertainty in A and Ai due to unknown and agA/mft 

contributions to the moments (see text for description). For the purpose of illustration, perfect 
experimental measurements of {sh — i^'d) = 0.30 GeV^ and {{sh ~^'d)'^) = 0.96 GeV^ have been 
assumed. The solid and dashed lines show the constraints on A and Ai from the first and second 
moments, respectively, while the area in the shaded ellipse shows the estimated allowed range. 
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FIG. 5. Estimate of the theoretical uncertainty in A and Ai due to unknown contributions 
from the shape of the electron spectrum, using the same approach as in the previous plot. The 
central value at 0(l/m^) is given by the black dot, while the dashed ellipse is the plot from the 
previous figure, shown here for comparison. Only the relative sizes and orientations of the two 
ellipses are meaningful in this figure, not their relative positions. 



13 



